Dynamics of Quintessence with Thermal Interactions 
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The cosmological dynamics of minimally coupled scalar field that couple to the background matter 
with thermal interactions is investigated in exponential potential. The conditions for the existence 
and stability of various critical points as well as their cosmological implications are obtained. Al- 
though we show that the effects of thermal interaction such as depressing the equation-of-state 
parameter of quintessence, is only important at the early time, the evolution of equation-of-state 
parameter of quintessence is manifested. The upper bound is required on the coupling between 
quintessence and relativistic relic particles such as photons and neutrinos. 
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Astronomical measurements from Supernovae la 0, 
galaxies clustering, for example, Sloan Digital Sky Sur- 
vey (SDSS) and anisotropies of Cosmic Microwave 
Background Radiation (CMBR) Q independently sug- 
gest that a large fraction of the energy density of our 
Universe is so-called dark energy, which has negative 
pressure, or equation of state with w = p/p < 0, and 
accelerates the expansion of the Universe. The origin of 
the dark energy remains elusive from the point of view 
of general relativity and standard particle physics. One 
candidate source of this missing energy component is a 
slowly evolving and spatially homogeneous scalar field, 
referred to as "quintessence" with w > -1 and 
"phantom" with w < — 1 0, H, El E| respectively. Since 
current observational constraint on the equation of state 
of the dark energy lies in the range —1.38 < w < —0.82 
[Tl| , it is still too early to rule out any of the above candi- 
dates. To study the global property of the cosmological 
system containing dark energy, phase space analysis is 
proved to be a powerful tool. 

Because dark energy redshifts more slowly than ordi- 
nary matter or radiation, there exists an important prob- 
lem that various proposals of the dark energy should ex- 
plain, that is, why the energy density of matter and dark 
energy should be comparable at the present epoch. A 
form of quintessence called "tracker fields", whose evo- 
lution is largely insensitive to initial conditions and at 
late times begin to dominate the universe with a nega- 
tive equation of state, was introduced to avoids the above 
problem |l2j| . Another approach to solve the puzzle is 
introducing an interaction term in the equations of mo- 
tion, which describes the energy flow between the dark 
energy and the rest matter (mainly the dark matter) in 
the universe. It is found that, with the help of a suitable 
coupling, it is possible to reproduce any scaling solutions. 

The effects of thermal coupling between the 



quintessence field and the ordinary matter particles 
has been investigated l.'i . Hsu and Murray set the 
quintessence field to be a static external source for a 
Euclidean path integral depicting the thermal degree of 
freedom and let the time-like boundary conditions of 
the path integral have a period which is decided by the 
inverse of temperature. They show that if in the early 
universe matter particles are in thermal equilibrium, 
quantum gravity [l4| will induce an effective thermal 
mass term for quintessence field </>, which takes the form 
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where T = $/p r o/ a is the temperature of the universe, 
Mp the Planck mass scale, (3 a dimensionless constant, 
a the scale factor (with current value ao = 1) and p r o 
the energy density of radiation at the present epoch, and 
find that even Planck-suppressed interactions between 
matter and the quintessence field can alter its evolution 
qualitatively. For convenience, Hsu and Murray |13| use 
an approximation that the back reaction on matter of 
quintessence is neglectable. However, this back-reaction 
effect, in fact, is significant both at the early time and late 
time. In this letter, we investigate the dynamics of the 
cosmology with quintessence using the phase space anal- 
ysis, which has the above thermal coupling to the matter 
in a complete manner. We shall seek the conditions for 
the existence and stability of various critical points as 
well as their cosmological implications. As will be shown 
in this paper, although the effects of thermal interac- 
tion, such as depressing the equation-of-state parameter 
of quintessence, is only important at the early time, the 
evolution of equation-of-state parameter of quintessence 
is manifested. The upper bound is required on the cou- 
pling between quintessence and relativistic relic particles 
such as photons and neutrinos. 

We first study the phase space of quintessence with a 
thermal coupling to ordinary matter particles in spatially 
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flat FRW cosmological background 



dt 2 -a 2 (t)dx. 2 
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For the spatially homogeneous scalar field minimally cou- 
pled to gravity with thermal interaction Q , the evolution 
is governed by the Klein-Gordon equation 



4> + 3H4> + V'(<f>) 
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where the overdots denote the derivative with respect to 
cosmic time and the prime denotes the derivative with 
respect to <f>. Here the Hubble parameter H 
determined by the Friedmann equation 
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where k 2 = 8ir/M 2 , p m and p m are the energy den- 
sity and pressure of the barotropic matter. According 
to Eqs.Q-lJSJ an d the conservation of energy, p m satis- 
fies the following continuous equation 



3H(p m+Pm ) = -2 (JL 
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and p m = (7 — l)pm, where 7 is a constant, < 7 < 2, 
such as radiation(7 = 4/3) or dust (7 = 1). It is clear 
that when the thermal coupling parameter j3 becomes 
zero, the equations ©-© will return to those of the 
standard one scalar field quintessence scenario. In cur- 
rent situation, both quintessence and barotropic matter 
are not conserved, but the overall energy is conserved. 
The energy density ratio of quintessence and barotropic 
matter satisfies the following conditions 

= r ( Q- + 9. _ 3 # 7 + 3H 14 ) , (7) 

\Pm P4> ) 

where the parameter Q = 2((3 / 'M p ) 2 (j)(j)p 1 ,o / a 4 , 70 = 1 + 
Wtfn, in which is the so-called parameter of state of 
quintessence which is defined by w$ = p^/pj,, and 
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Rewriting the equations © and ©, we have 
Pm + 'SH(w e n { f + l)p m = 0, 
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where wffi =7 — 1 
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and tUj, = 



3M£Ha 4 pr, 

ZM'^Hatp ' respectively. Now, introducing the following 
variables: 
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the equation system©-© becomes the following system: 
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The energy density parameter and the equation-of-state 
parameter of quintessence, fi^ and w^, satisfy the con- 
straint equation 
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and — (x 2 — y 2 )/(x 2 + y 2 ), respectively. Since the 
energy density of the barotropic fluid p m is semi-positive- 
definite, any cosmological model can be represented as a 
trajectory in the phase space that is bounded within the 
unit disc, i.e. < x 2 + y 2 < 1. From the definitions 
of these new variables, It is not difficult to reduce the 
effective parameter of state of quintessence to 
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To be concrete, we consider the quintessence with an 
exponential potential energy density, i.e., 



V(4>) = V exp(-A K(? 



(22) 
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where the parameter Vq and Ao are two constants. Expo- 
nential potentials have been studied extensively in var- 
ious situations, and these are of interest for two main 
reasons. Firstly, they can be derived from a good can- 
didate of fundamental theory for such being string/M- 
theory; Secondly, the motion equations can be written 
as an autonomous system in the situation. The au- 
thor of an earlier work considered a scalar field 
with a single exponential potential in a homogeneous and 
isotropic universe. Subsequently, it was extended to in- 
clude barotropic matter [Ty,[l3 and multiple scalar fields 
|l8j. and generalized for anisotropic universe Ac- 
cording to the definitions, the parameters A and T both 
become constants and are equal to Ao and 1, respectively. 
Under the circumstance, the equations l|19|l constitute an 
autonomous system as follows, 

% = \y[l{l~x 2 ~y 2 ) + 2x 2 ]~^\ a xy, 
dz 3 
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From the equation of state of matter © and Eqs. I|12|l- 
IjlSp. we obtain that 
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where we have chosen the exponential potential (|22l) . p 71 
can be directly expressed as 
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where the parameter pm,o is the present value of energy 
density of matter. It is easy to find that in the case of 
/3 = 0, the energy density of matter will evolve in the 



manner of p n 



for 7 = 1 as usual. 



When barotropic matter is under consideration and/or 
the equations of motion are too difficult to solve analyt- 
ically, phase space methods become particularly useful, 
because numerical solutions with random initial condi- 
tions usually do not expose all the interesting properties. 
In table U we list the critical points, conditions for their 
existence and the cosmological parameters there. 

To investigate the stability of these critical points, we 
can write the variables near these points (x c ,y c ,z c ,£c) m 
the form x — x c + rjx,y = y c + V2, z — z c + 773 and 
£ = £ c + r)4 with rji,i = 1, 2, 3, 4, the perturbations of the 
variables about the critical points to the first order. This 
leads to the equation of motion, 



U' = A U 

where the 4-column vector U = (r]i) T , 



(26) 



i = 1, 2, 3, 4 rep- 
resents the perturbations of the variables and A is a con- 
stant 4x4 matrix. For stability we require the all 4 
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FIG. 1: The portrait of scaling behavior of the quintessence 
with thermal interaction for exponential models V(4>) = 
Vo exp( — AqK0), where Ao = 25 and the parameter 7=1. 



eigenvalues of A to be negative. For the critical points 
listed in table 1, we find the eigenvalues of the linear 
perturbation matrix for different stable critical points, 
see table ITT1 

Critical points of type A correspond to the scalar field 
dominated solutions (O^, = 1), which exist for sufficiently 
flat potentials, A 2 < 24. Moreover, for the barotropic in- 
dex 7 > A 2 /18, they are stable to the extent that this 
kind of solutions represent late-time attractors in the 
presence of a barotropic fluid. Note that type E is a 
special case of type A, i.e., A = 0. 

Critical points of type B correspond to another kind of 
late-time attractor where neither the quintessence field 
nor the barotropic fluid completely dominates the evo- 
lution of the universe. They are known as scaling so- 
lutions where the energy density of quintessence is pro- 
portionated to that of barotropic fluid at late time with 
17^ = I87/ A 2 . It is remarkable that the barotropic in- 
dex 7 is required to be no more than 4/3 for the sake of 
stability. 

Critical points of type C denote the barotropic fluid 
dominated solutions where 17^ = 0. They are unstable 
for all reasonable values of 7 and A (i.e. however steep 
the potential). The ones of type D are the so-called 
kinetic-dominated solutions where the late-time evolu- 
tion of the universe is dominated by the kinetic energy 
of quintessence with a stiff equation of state, i.e. = 1. 
These solutions are unstable as one expects. 

As shown in table [I] all the critical points exist when 
and only when z — 0, that is to say that f3 does not ap- 
pear in the final expressions that determine the critical 
points. This means that with the evolution of the uni- 
verse the thermal coupling between the quintessence 
and matter always goes asymptotically to vanish away. 
Fig^ shows a scaling solution that the quintessence goes 
towards an attractor. 

In the above, we have studied the phase space of scalar 
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Type 


critical points (x c ,yc,z c ,^ c ) 


A 






stability 


A 


i ±\A - & o. 


any 


1 


-1 + T8 


stable for A 2 < 24 and 7 > 4§~ 


B 

C 
D 

E 


3 J, ±3^^1,0,0 
0,0,0,any 
±1,0,0,0 
0,±1,0,0 


A#0 
any 
any 

A = 


187 

as 



1 

1 


7-1 
undefined 
1 

-1 


stable for 7 < | and A 2 > I87 

unstable 
unstable 
stable 



TABLE I: The critical points and their physical properties for models with exponential potential. 
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Critical points (x c , Vc, z c , £ c 



eigenvalues 



A 
B 



C 
D 

E 



±x 1 



f ,0,0 



32, ±3^1^,0,0 



0,0,0, any 
±1,0,0,0 
0,±1,0,0 
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2)- J (2-7) (2-97+^) 



(7 - 2) + J(2- 7 ) (2-97 
0, 



± 
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Q I 37 37 

2 ' 2 ' 



-2 ± ^ 

0,6 -3 7 , 1,3 T I 2 



0,-3, 



-3 7 



TABLE II: The eigenvalues of the critical points for the exponential potential. 



field with thermal interactions to cosmic matter in a flat 
FRW cosmological background. The critical points in- 
dicate that these stable attractor phases corresponding 
to the vanishing of the thermal interaction between the 
quintessence and the barotropic matter. However, what 
will be the effects of the thermal interactions in the cir- 
cumstances? In fact, the introduction of the interacting 
term will change the evolutionary tracks of the dynami- 
cal systems, which was not well manifested by the above 
qualitative analysis. Therefore, we study their dynamical 
evolution numerically and compare the results with the 
cases of no coupling quintessence scenario. In Fig|21 the 
evolution of equation-of-state parameter of quintessence 
is manifested. Obviously, the effective equation-of-state 
parameter of quintessence vfl that incorporate the ef- 
fect of the thermal coupling between quintessence and 
matter is smaller than that of quintessence without cou- 
pling to the background matter. 

In conclusion, we have investigated the cosmological 
dynamics of scalar field, quintessence, that couple to the 
background matter with thermal interactions in exponen- 
tial potential. The conditions for the existence and stabil- 
ity of various critical points are obtained by phase space 
analysis. We find that incorporating the effect of ther- 
mal coupling between the quintessence and background 
matter does not qualitatively alter the late-time evolu- 
tion of the components in the universe. Just as the or- 
dinary situation that there is no direct coupling between 
quintessence and background barotropic matter ^|| > f° r 
the parameters A 2 < I87, the universe is dominated by 
quintessence at late time; and for A 2 > I87, quintessence 
does not entirely dominate the universe but remains a 
fixed fraction of the total matter at late time. However, 
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FIG. 2: Evolution of the effective parameter of state w^f for 
j3 = 1.5, /3 = 1 and j3 = 0, respectively, where we let 7 = 1, 
Aq = 2 and z in the label of horizontal axis denote red shift. 



we also find that the effects of thermal interaction make 
the equation-of-state parameter of quintessence become 
lower than the one in the standard quintessence scenario 
of dark energy at the early time and depress the energy 
density of matter simultaneously, see F'igO I n other 
words, large coupling to relic particles such as neutrinos 
can be ruled out as they lead to a problematic equation- 
of-state parameter, which is consistent with Ref. |13j . 
The fact that <j> is probably close to the Planck energy 
suggests that more interactions should figure in to the be- 
havior of quintessence. However, this will possibly lead 
to a relic density problem at nucleosynthesis. Finally, it 
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is worth emphasizing that this thermal interaction may 
be helpful to understand the coincidence of fl m ~ 0^. 
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FIG. 3: Evolution of the energy density of barotropic matter 
p m for f3 = 1 and /? = 0, respectively, where we let p m _o = 1, 
7 = 1, A = 2 and a in the label of horizontal axis denote red 
shift. 



This work is supported by Shanghai Municipal Ed- 
ucation Commission No. 04DC28, Shanghai Munici- 
pal Science and Technology Commission No. 04dz05905 
and National Natural Science Foundation of China under 
Grant No. 10473007. 



[1] A. G. Riess et al, Astron. J. 116 (1998) 1009; S. Perl- 
mutter et al, Astrophys. J. 517 (1999) 565; J. L. Tonry 
et al, Astrophys. J. 594 (2003) 1. 

[2] M. Tegmark, et al, Astrophy. J. 606 (2004) 702. 

[3] D. N. Spergel et al, Astrophys. J. Suppl. 148 (2003) 175. 

[4] P. J. E. Peebles, B. Ratra, Rev. Mod. Phys. 75 (2003) 
599; T. Padmanabhan, Phys. Rep. 380 (2003) 235. 

[5] R. R. Caldwell , R. Dave and P. J. Steinhardt, Phys. Rev. 
Lett. 80 (1998) 1582. 

[6] X. Z. Li, J. G. Hao and D. J. Liu, Class. Quantum Grav. 
19 (2002) 6049. 

[7] R. R. Caldwell, Phys. Lett. B545 (2002) 23; R. R. Cald- 
well, M. Kamionkowski and N. N. Weinberg, Phys. Rev. 
Lett. 91 (2003) 071301. 

[8] S. M. Carroll, M. Hoffman and M. Trodden, Phys. Rev. 
D68 (2003) 023509. 

[9] J. G. Hao and X. Z. Li, Phys. Rev. D67 (2003) 107303; 
Phys. Rev. D68 (2003) 043501; Phys. Rev. D70 (2004) 
043509; D. J. Liu and X. Z. Li, Phys. Rev. D68 (2003) 
067301. 

[10] P. Singh, M. Sami and N. Dadhich Phys. Rev. D68 
(2003) 023522; S. Nojiri and S. D. Odintsov, Phys. Lett. 
B562 (2003) 147; B571 (2003) 1; P. F. Gonzales-Diaz, 
Phys. Rev. D68 (2003) 021303; L. P. Chimento and R. 
Lazkoz, Phys. Rev. Lett. 91 (2003) 211301; H. Stefancic, 



Phys. Lett. B580 (2004) 5. 
[11] A. Melchiorri, L. Mersini, C. J. Odman and M. Trodden, 

Phys. Rev. D68 (2003) 043509. 
[12] I. Zlatev, L. Wang and P. J. Steinhardt, Phys. Rev. Lett. 

82 (1999) 896; P. J. Steinhardt, L. Wang and I. Zlatev, 

Phys. Rev. D59 (1999) 123504. 
[13] S. Hsu and B. Murray, Phys. Lett. B595 (2004) 16. 
[14] M. Kamionkowski and J. March-Russell, Phys. Lett. 

B282 (1992) 137; R. Holman, S. Hsu, T. W. Kephart, E. 

W. Kolb, R. Watkins and J. March-Rusell, Phys. Lett. 

B282 (1992) 132; S. M. Carroll, Phys. Rev. Lett. 81 

(1998) 3067. 

[15] J. J. Halliwell, Phys. Lett. B185 (1987) 341. 

[16] E. J. Copeland, A. R. Liddle and D. Wands, Phys. Rev. 

D57 (1998) 4686. 
[17] R. J. van den Hoogcn, A. A. Coley and D. Wands, Class. 

Quantum Grav. 16 (1999) 1843; I. P. C. Heard and D. 

Wands, Class. Quantum Grav. 19 (2002) 5435. 
[18] X. Z. Li and J. G. Hao, Phys. Rev. D69 (2004) 107303; 

J. G. Hao and X. Z. Li, Class. Quantum Grav. 21 (2004) 

4771; Phys. Lett. B606 (2005) 7. 
[19] A. P. Billyard, A. A. Coley, R. J. van de Hoogen, Phys. 

Rev. D58 (1998) 123501. 



